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1. Introduction 


This data report is a supplemental record of data and calcu- 
lations related to the M.S. thesis of Douglas L. Huizenga, 
Protonation Characteristics of the Dissolved Organic Matter in 
Seawater (University of Rhode Island, 1977). The individual 
fraction titration data, calculated organic matter titration 
curves, model parameters, and copies of computer programs are 


presented. 


Organic Matter Samples 


The organic matter samples investigated are listed in Table 1. 
The location, date, and depth of sample collection are given. 
A specific blank is associated with each sample and is given in the 
last column of Table 1. Blanks were prepared to be equivalent on 
a volume basis with the samples. The last ten samples were all 


handled in the same manner and only have a single blank sample. 


Titration Names 


Each titration was assigned a unique name, such as GGI371A. 
The first letter(s) plus the first number ("GG1") indicates the 
sample code and fraction number (here sample GG, fraction 1). 
The rest of the name specifies the lab book page number for the 
titration (here page 371A). A few titrations of 0.7m NaCl0, 
without added sample were performed. The names for titration of 


medium alone begin with "MD", A fraction number of 0 indicates 


that the sample was not obtained by the charcoal isolation procedure. 
A fraction number of 4 indicates a sample where all three fractions 


from charcoal were combined. 


Titrations Performed 


Titration information for every titration is given in Table 2. 
The order of samples in this table is also the order in which the 
titration data are presented. The use of the titration parameters 
is given in Appendix A of the thesis. All the information about 


a titration that is needed to do the calculations is given in this 


table, 
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Titration information. 


Type, OM is an organic matter sample 
BL is a blank sample 
w, weight of water in sample at start of titration, g 


[OM], amount of organic carbon in sample at start of 
titration, mg 


E , pH electrode calibration E , mV 


0 grams of titrant solution / ml 


t’ 
“H0? grams of H,0 / gram titrant solution 

OH,» moles of OH / kg titrant 

Date, date when titration was performed 

Blank, blank appropriate for that sample 

Sample Weight, weight of sample taken initially, g 


I, ionic strength of solution at start of titration, m 


Factor, that fraction of initial sample in the test solution 
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II. Acid-Base Titration Data 


bi i a a a te Ra dl et slit) 


The data for titration of organic matter samples and blanks 


are presented in Table 3. 
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Table 3. Experimental data from titrations. Ordering of data 


is nearly that of Table 2. 


ML, volume of titrant added, ml 


EH, potential of electrode pair, mV 
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38. 


III. Organic Matter Titration Curves 


Calculations for each titration made by the program TCRVDI 
are reported in Table 4. Weighted averages of the seperate fraction 
titration curves were made, producing combined titration curves. 


These combined titration curves are reported in Table 5S. 


The averages were made with values interpolated to each 0.1 pH unit. 
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Calculated organic matter titration curves for organic 


sample - blank pairs. Ordering as in Table 2. 


PH, calculated pH 


MCMPOC, OH consumed by organic matter, uwmoles/mg OC 
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Combined titration curve data. 


Units are umoles/mg OC. 
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IV. Model Calculations for the Type 1 Sites 


The modified Henderson-Hasselbach equation was used as the basis 
for non-linear least squares fits of the Type 1 site titration curves. 
Calculations were performed over the pH range 3 to 6. Model parameters 
for fits calculated by program NLNOM for the seperate fraction titration 
curves are given in Table 6. Calculation results with NLNOM for the 
combined titration curves are given in the thesis. Calculations 
performed subsequent to the thesis with the SAS 76 NLIN program for 
the combined titration curves are given in Table 7. 

The modified Henderson~Hasselbach equation as used for the 


calculations was: 


X = €2/(10**( (pK, -pH) /N) +1) -C1 


where C2, total number of sites or T(-COOH) 
a constant (that pH where half the sites are dissociated) 
negative log of the hydrogen ion concentration 
a constant 


MCMPOC or the umoles of base consumed/mg OC 


Table 6. Parameters from model fit to Type 1 sites for the 
individual fraction titrations calculated with 


) program NLNOM. 
4 


N, a constant 
PK,> a constant (that pH where half the sites are dissociated) 
| T(-COOH), C2 or the total number of Type | sites 


if Cl, adjustable offset of the titration curve 


SUMSQ, sum of squares, £(X 


| 
| 
if 
] 


) 


n, number of data points used in the calculation 


observed “estimate 


2S gee eas Do 


73 


| 
{ 
1 T(-COOH) 
| Titration N PK, umoles/mg OC Cl SUMSQ n 
: LL CLL j 
| GG1371A 1.99 3.59 13.85 3.28 0.0413 16 
| GG1392A 2.10 3.38 14.82 2.91 0.0224 16 
1} GG2373A 1.64 3.92 4.10 1.36 0.0202 15 
{ GG2421B 1.60 4.10 3.39 1.67 0.0223 10 
| GG3374B 1.60 3.73 6.09 -1.05 0.0253 16 
GG3421B 1.70 3.80 5.08 0.70 0.00275 15 
| HH1377B 1.91 3.66 14.86 0.93 0.0356 16 
' HH1399C 2.32 3.13 18.41 4.32 0.0433 16 
HH2378B 1.86 3.90 4.56 0.48 0.0216 14 
i HH2399B 1.66 4.04 3.65 0.67 0.0261 13 
: HH3378B 1.88 3.67 6.88 0.20 0.0176 16 
4 HH3399B 2.08 3.61 7.05 1.66 0.0087 13 
4 111380A 1.74 3.70 12.50 0.37 0.0599 17 
‘ 111399C 1.76 3.76 11.57 1.04 0.0321 15 
1 112381B 1.41 4.14 3.70 0.10 0.0115 15 
| 112404a 1.39 4,08 3.41 0.83 0.0072 12 
: 113382B 1.81 3.90 7.01 1.27 0.0832 16 
i 1I3404A 2.29 3.46 9.20 1.18 0.0044 12 
JJ1383B 1.11 4.68 22.28 -3.37 1.041 16 
JJ1407A 1.07 4.68 21.23 -0.66 0.852 16 
JJ2386A 1.74 3.94 4.47 0.04 0.0096 15 
JJ2405A 1.43 4.07 3.83 -0.31 0.0032 13 
JJ33874 2.02 3.51 8.15 1.32 0.0136 15 
JJ3405A 1.80 3.71 7.14 -0.11 0.0127 9 
KK1383B 2.53 3.39 15.95 2.52 0.0388 12 
KK1407A 2.94 3.15 17.88 5.88 0.0357 15 
KK2386A 1.40 4.12 3.79 0.59 0.0018 14 
KK2408A 1.36 4.18 3.58 2.46 0.0176 14 
KK3387A 1.97 3.72 7.74 -3.32 0.0129 15 
KK3408A 2.07 3.83 7.11 2.16 0.0147 14 
LL1392A 1.28 4.06 5.03 -0.68 0.0091 13 
LL1410A 1.32 3.96 5.33 0.04 0.0250 15 
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T(-COOH) 

Titration N PK, umoles/mg OC Cl SUMSQ n 
LL2393A 182 4.85 7.06 0.92 0.0394 ‘7 
LL2409A 1.70 4.79 6.56 2,22 0.0315 13 
LL3393A 1.92 3.94 6.26 0.77 0.0049 lo 
LL3409A 2.20 3.60 7.89 2.04 0.0096 12 
MM1396A 2.85 4.34 15.66 3.40 0.0264 17 
MM1410A 2.73 3.96 13.61 4.05 0.0364 17 
MM2394B 1.84 3.86 5.09 0.30 0.0096 15 
MM2411A 1.54 4.07 3.94 1.56 0.0049 14 
MM3394A 1.70 3.80 6.43 1.43 0.0085 14 
MM3411A 1.94 3.83 7.59 2.09 0.0147 14 
NN1396A 2.28 3.42 13.29 1.49 0.0090 1S 
NN1413A 2.01 3.65 11.13 1.75 0.0301 16 
NN239SA 1.62 4.00 4.04 1.16 0.0026 15 
NN2412A t.S2 4.10 3.86 1.22 0.0062 1S 
NN339SA 1.87 See.) 5.32 1.32 0.0037 14 
NN3412A 1.62 3.83 $3.73 -0.66 0.0159 1S 
001413A 2.09 3.85 13.96 1.85 0.0135 15 
001423A 2.09 3.54 14.36 1.95 0.0081 uF 
002414A 1.79 3.95 6.45 0.71 0.0029 13 
002420A 1.82 3.93 6.12 1.26 0.0026 17 
003414A 1.94 3.66 8.42 1.54 0.0077 lo 
003420A 1.95 3.66 8.39 1.30 0.0028 1S 
QQ1423A 2.46 3.10 17.43 5.74 0.0138 13 
QQ1426A 1.93 3.83 12.83 2.78 0.0142 16 
QQ2424A 2.23 3.54 6.94 2.52 0.0055 12 
QQ2425A 1.73 3.96 5.17 1.38 0.0085 14 
QQ3424A 1.89 3.71 7.44 2.09 0.0052 14 
QQ3425A 2.28 3.35 9.65 3.43 0.0220 13 

10. SHrUVQQ1433-2 1.58 3.92 21.65 3.40 0.1724 15 
SHTUVQQ1433-4 1.77 3.61 18.30 2.17 0.0243 13 
LHrUVQQ1433-6 1.92 3.50 15.39 2.91 0.0104 11 


T(-COOH) 
Titration - umoles/mg OC Cl 


N1237A : : 1S. -99 
N1255B 3. " 22. .79 
N22S0B ; ‘ 6. .00 
N3247B 

P1318A 

P1340B 

P2320A 

P2338 


P3322 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 


no sil 


Table 
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Parameter values and asymptotically valid errors for 
selected combined titration curves calculated with 


SAS 76 program NLIN. 


TCOOH, C2 or the total number of Type 1 sites 


PK, a constant (that pH where half the sites are dissociated) 


N, a constant 


C, adjustable offset of the titration curve 


YY 
_— <5 wee Estimate Std. Error 95% Confidence Interval 
Sample N 
TCOOH 13,815 0,267 13.268 14.363 
pk 3.506 0,026 3.454 3.559 
N 1.851 0,031 1,786 1.916 
(2 4,781 0,237 4,294 5.268 
Sample p 
TCOOH 15,509 0.276 12.943 14.075 
pk 3,328 0,029 3,268 3.3588 
a 
N 1.955 0,030 1.892 2.018 
| C 5.434 0.252 4.917 5.950 
Sample 6G 
TCOOH 9,730 OLehi 9,297 10.164 
| 
pk 3.606 0.028 3.548 3.664 
N 1,924 0,039 1.8435 2.008 
C 3.146 0.182 2.773 3.520 
Sample HH 
TCOOH 10,89] 0.31] 10,253 11.530 
pK 3.575 0.038 3,496 3.653 
N 1.982 0,052 1.874 2.089 
C 5,661 0,269 3.109 4.212 
a 11 
wad 0,596 0.247 9,088 10.103 
"47 0,030 3.685 5,810 
a 0.0468 1.680 1,878 


or 2.168 3.017 


brit Sse pe! 


ln 


Estimate 


8,963 
3,666 
1,982 
2.809 


10.331 
3.674 
1,965 
3.214 


10,025 
3.557 
1.982 
3,449 


Std, Error 


0.276 
0.035 
0,052 
0.231 


0.121 
0.017 
0.026 
0.102 


0,134 
0.016 
0.025 
0.113 


0.163 
0,022 
0.029 
0.142 


7B 


95% Confidence Interval 


11,049 
3.491 
2.273 
3.782 


8.714 
3.631 
1,928 
2.600 


10.056 
~ 3..640 
1,923 
2.983 


9.690 
3.512 
1,922 
3.158 


12.184 
3.634 
2,487 
4.732 


10.606 
3.708 
2.017 
3.446 


10.360 
3.602 
2.043 
3.740 


| 
| 
| 
| 


da inte IS 


aad 
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V. Programs Used for Calculations 


TCRVDI and NLNOM 

These programs were used on the laboratory Varian computer 
system. TCRVDI performs curve difference calculations and was used 
to subtract blank titration curves from organic sample titration 
curves. NLNOM performs the least squares fit of the model to the 


titration curve over a specified pH range. 


SAS 76 NLIN 


A copy of the cards used for running the SAS 76 NLIN program 


are given. The Marquardt method was used for all calculations. 
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a eh ee 


PAGE 1 08/@9/77 VORTEX CTUF 826 HOURS | 
1 TITLE TCRVDI 
i 2 CHHHHAHREHHHHAEHAEREARHHEHRAAHENHRHEHHER RAHA AR RRR 
3 ¢ THIS MODULE WILL CRUNCH 
i 4 ¢ TEST CRVDIF DATA ONLY IT IS PROTO 
{ 
6 6 
| Z  € THIS IS REVISION B : THIS REVISION. INCORPORATES FLEXIBLE. 1/0 
| $ ¢ ~- ALSO IT SUPRESSES THE IMAGE DUMP, 
| 9 ¢ AND RATHER OUTPUTS ONLY THE VITAL DATA | 
| 1 ** ALSO NOTE THAT THE VRSION OF | CHVDIF | 
| ae USED HERE IN IS THE MOST CURRENT REVISION | 
| i2 ¢ AND SHOULD BE THE VERSION USED IN THE GLOBAL : 
| 13 —¢ OMTITR SYSTEM. SR | 
| mu ¢ ### REV B INCORPORATES SIGN OFF FROM ROOT SEOWENT | 
1s ¢ KICKS OUT LAST P 
17 ¢ REVISION C: | THIS REVISION CORRECTLY INTERPOLATES | 
ig ¢ DATA SETS WHICH DO NOT HAVE MATCHING ENDPOINTS 
19 «C¢ (IF NO INTERPOLATION IS POSSIBLE THE PROGRAM 
2 2 ¢ SO NOTIFIES THE OPERATOR. ADDITIONALLY THE 
21 ¢ PROGRAM NOW SIGNS’OFF TO INDICATE COMPLETION. 
| 
¢ RIFFLE 28SP76 (SEE OMTITR FOR SYSTEW 
: C 
¢ 


22 
23 
24 
a 

27 C CoeeoosoooooseesooeoeeeNOTE EXTERNAL STMTS ARE RESTRICTED TO ROOTONLY 

28 EXTERNAL BPFCBS 

29 COMMON /FCBS/ IBFCB (13), IBLU, IPFCB (13), IPLU /INTERA/INTER 

30 COMMON /IOSPEC/ LISTLU, IOLUN 

31 Coecccccccecccncccocevencocese NOTE ALL ROOTS MUST MAKE ONE WRITE STMT 

32 Cc THIS PUTS VSFORTIO IN PROG, A MUST 
( 33 . —Laectiicek  €6€6ClCCC:C:*;~”:”” 

34 C __LET HIM KNOW WHAT’S GOING O 

35 (oseseccsvenncecosconnassocevennenssooeennnasooeseaseoooesnssoooees 
36 
37 
38 
39 
40 
41 
42 
43 
44 


= 


Cc SPECIFY PROGS NAME AS PER TITLE 


WRITE( 1,920) 
920 FORMAT(’ {THIS IS *°eTCRVDI REV cee’) 
oo ; RESPECFICATION OF ted 
M CURRENT 


Cc IOLUN FRO 
Cc (NOTE THIS IS ALWAY FIRST DONE AT OC. 
45 GO TO 921 
46 8040 WRITE(IOLUN ,922) 
47 READ(IOLUN,923) LISTLU, IOLUN 
pe c GO TO 924 


50 Cc 
$1 921 WRITE(!,922) 
$2 922 FORMAT(? SPECIFY pun IOLUN'/' +### +###') 


$3 READ (1,923) LISTLU, I0L 
$4 923 FORMAT (1X, 14, 1X, 14) 
SS C..+e ree 

S56 C ECHO 

57 


Coss cssccecs 
58 924 WRITECIOLUN ,920) 


in MN ii > ik) = Aa 


ee 


2 
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88/09/77 VORTEX CTUF 824 HOURS 
2 WRITECIOLUN ,923) LISTLU, IOLUN 
C 
Cc 
Caeccwestrseseeccs stat eS Stee SSeS PES TESST SEC ESESESL SUK ESER ESS 
C PROG REALLY STARTS HERE 


Ceseteee est et ers eee Sr Ss SS SSS SE SESS SS SSERSESERESS RAE SES SERSEE 


$4 WRITE(IOLUN, 205) 

20S FORMAT (° vara TAPE ON DRIVE ?YES=1'/ 
+ % +g? 
READ (IOLUN , 206) IYES 

206 FORMAT (1X, 12) 
IF(IYES.NE.1) GO TO 84 


lolole) 


CALL CRVDIF 
WRITECIOLUN , 436) 
436 A aaa * s:::TCRVDI BYE®:::::°) 


Dd ddd ddd add ddd did dad dy dade A dad dudsds dad dudadsdidud ed dudadsdudud dadudad dadaddadad dudes dada dad cdadsd dada 
CRVDIF DOUG’S "CURVE DIFFERENCE" PROG 


THIS ROUTINE CALCULATES DIFFERENTIAL TITER UPTAKE 
OF SEA WATER SAMPLES BY ORGANIC MATTER, BASED ON 
POTENTIOMETRIC TITRATION WITH OH’ RE: WRITEUP 


THIS IS REVISION C: WHICH ENCORPORATES THE INTERACTIVE 
LOGICAL UNIT SPECIFICATION IN RUNTINE. INITIALLY 
INTEROGATION IS TO THE OPCOM DEVICE. 


THE Sniper eniee HAS CORRECTED TO HANDLE 
NONCOINCIDENT ENDPOINTS, BOTH ENDS, AND Here — 
IS NOTIFIED IF INTERPOLATION IS NOT POSSIB 


aaNAANANAAAANAANANAAAAA SS .gaANMAANAANMOAANANANANAN 


ALSO I HAVE OMMITED THE SUPERFLUOUS CORE IMAGE DUMP 
WHICH WAS NEEDED DURING DEBUG SEGMENT. SR. 


| 
| 


Mea Dy RA. Sue a RE LH Na Bao 
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PAGE 3 08/09/77 VORTEX CTUF $24 HOURS 
117, 
1s oC THIS PROG ASSUMES DATA BAS BEEN PLACED ON LUN 174=MT IN 
| lig ¢ THE FORMAT SET FORTH IN WRITE UPFOR CRVIF 
121 CG TO BE PROTOTYPE FOR 
122 € SUBROUTINE (OVLAY) 
124 : THIS PROG REOQUIRS CALLER TO INITILIZ.::: | 
1255 ¢ : COMMON BLOCK /FCBS/ VIA EXTLURNAL STMT:: : | 
136 C  REPRENCNG BLOCK ATA PROG BPECBS ei 
| 127 ¢ THAT PROG IS OM RESIDNT vA 
129 «¢ (CALLING PROG MUST HAVE AT LEAST ONE WRITE STWT)) 
130 (CIF THE PROG IS TO BE RUN AS OV)) 
1 ¢ 28SP76 RIFFLE 
133 SUBROUTINE CRVDIF 
1 134 C SET STORGE 
135 COMMON /IOSPEC/ LISTLU, IOLUN 
4 136 COMMON /FCBS/ IBFCB (13), IBLU, IPFCB (13), IPLU /INTERA/INTER 
; 137 DIMENSION RPLOT (15), IPLOT (66), LABLX (15), LABLY (15) 
j 138 DIMENSION DATA (60,6), NCODE (8)) NOTE (35), PARAM (12) 
: 139 DIMENSION MYFCB (13) 
140 DATA (RPLOT (1) ,1=1,9)/1.0,1.0,0.0,7.0, 1.0,1.0,+90.,9.0, 
i a +1. 
: 142 DATA (LABLX (1), I=1,15) /7®2H ,2HPB, 7@2H / 
143 DATA (LABLY (1) ,1=1, 15) /* M*, "1G", "RO", "MO", °LE","S °, 
144 1 °OH’,* C*,*ON')*SU', ME’, D/?, MG,’ 0°,°C ¢/ 
EE eerie AIA ib eI 2c A A AR A Alc ae es 
146 OC GET THE LUNS 
14? 4 DATA AREA ON 
14s DISC VIA BPFCBS 
149 (MAIN INITIALIZS) 
ise 
151 LUN* IBLU 
152 BDO. 1 [#1,13 
| 153 MYFCB (1) =iBFCB (1) 
154 1 CONTINUE 
| 155 IREW=0 
156 CALL V8OPEN (13,LUN,MYFCB, IREW) 
is? C 
iss 
10 Gece enmceamnisnmniensinninainitltonasnctalindichtinildnetlinalid pine entnenrts terme 
ie Ci SET LOOP PARAMS 
161 ¢ WIL BE 2 PASS’S 
162 ISET=1 
163 ICLEAR=6 
164 NSET=1 
165 INDEXe | 
166 Go TO 3 
167 OC 
168 2 ISET=ISET+1 
169 ICLEAR=4 
170 NSET=5 
171 INDEX=7 
1 RI Ng hs aad Wh Fl. 0" TROVE 23. 5 .. FIRST TIM READ PLT HEDR 
173 3°" iF CSET BO! 1) READ’ (174, °4)  CiPLOT’ (i) ; 148,60) 
174 4  FORMAT(IX,2A2, 17,342,' T14,4A2, 12374 


a 


| 
| 


68/89/77 VORTEX 824 HOURS 


READ PARAMS EACH D SET 
READ(174,5) NCODE (NSET), NCODE (NSET+1), NCODE(NSET+2) , 
1 NCODE (NSET+3) 
FORMAT (1X, 4A2) 
READ (174, 62) (NOTE(I) , I1=1,35) 
60 FORMAT (1X, 3SA2) 
paciaae: 61) PARAM( INDEX) » PARAM (INDEX+ 1) , PARAM ( INDEX+2) , 
1 PARAM(INDEX+3) , PARAM ( INDEX+4) , PARAM ( INDEX+5) 
61 FORMAT(1X,F7.3, T10,F5S.2, T16, F7.2, 124,E11.4, 
1 1T36,E11.4, T48, E11.4) 


$1 J=1, ICLEAR 
DO 6 N=1,60 
DATA (N,J)=0.0 
CONT INUE 
CONT INUE 
READ RAW DATA 
ie ae us ) VOLUME OR 68 PTS 
HE CURRENT LIMITATION 
NPAIRS=0 
NPAIRS=NPAIRS+ 1 


DATA (NPAIRS, 1) =VOLML 
DATA (NPAIRS,2) =EH 
IP (CNPAIRS.EQ.60) .OR. (VOLML.LT.0.0)) GO TO 9 


KILL FLAG ROW 
(FLAG IS ANY -VOLML) 
( THE EH IS ARBITRARY&#0) 


COCSOSSHEHSSOHESHHOOHHEEHLOSOELOOSOHSEOHOOOOHOSOSHES HOSE SOE OOSOEOOOOEOOOE 


ORDER WRT ML (INC DWN) 
RDER3 OM RESIDNT 


is) 
COPSOSSOSSEESSSHE SESH OS OSEOSEOOSEOESHSOSOSOHOOOOEHOEEOOOOSEEOOOOOEOOS 


MXRSUB=NPAIRS 
MXCSUB=2 
NCOL=1 
11DC6=1 


tae ee | | | | | | ce | eC | ee | cr | cn | | eer ee em ym SY 
i i ee A a a oo a ee 


THIS SEGMENT IS PURPOSELY MASKED, FOR FUTURE TEST 


WRITE(LISTLU, 1848) ((DATA (J,K) ,K=1,6), J=1,NPAIRS) 
FORMAT(* °,T2,OPF7.3, 1T10,@PF7.2, T1T18,QPF8.4, T27, 1PE11.4, 
1 T39,OPF8.4, TS@, 1PE11.4) 


Uae a ee ee | eC | | ee a | a | | eT er | 
es es fs FS 8 2 8 8 8 8 See se Ss Ss se se se esse "s"s"e"e"s"se"e"e 


agNAANNAAAANNN 


i Ae Nia nt aha, 


84, 
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CE OE Ci tnth A ORR Oot RAN Be CAO OOOO AO AINE SET ROW COUNT IN CURRENT X 
IPLOT (10)=NPAIRS 


lolelelele) 


COCCOCOOOOEEOHOEESOOHEESOEEOOOOEEOOOCOOESOOOOO CES OOOO OS OOSEOSOOOSOOCOOS 
CALC MICROMO 


/ ML 
COSSOSOSSEOEHSSSOEHOSOSSSEOHHSHSSSOSSOSOSOSHOSOSOHESEOHSOSHOHEOHOOOOCEE 


pia a (INDEX+3)® PARAM (INDEX+5)®(10.¢*6) 


CONTI 
IF(ISET.EQ.1) GO T 
CONVRS= (PARAM (1) /PARAM (7) ) *CONVRS 


COOSOEEOEOSSOHEEOOEOOHSOOOCESOONSSOOEEEEOOOSSEOCOOOOEOOOOOOECOOSCE 
CALC PRELIM PH & CONSUMP 

COCCCOOCCOCEOOEECOE OEE OOEOOCOLOEHOOSOEESESOOOEESEOESEOHHEHOEOESOEDEECE 
76 DO 16 J=},NPAIRS 

DATA (J,3)= ~((DATA (J,2)-PARAM (INDEX+2))/59. 16) 

IF (ISET .EQ. 1) DATACS, 5) =DATA(J, 3) 

yee 4)=DATA(J, 1) 
CISET EQ. 1) DATA: 6) =DATA(J ,4) 


ABOQAA 


i 10 CONTIN 
3 c 
y ¢ 
| 
| ¢ STIK ROW COUNT IN UNUSED 
q ‘ IPLOT (38+ISET) =NPAIRS 

¢ 
| 
. COSSSSOOHOSOHSHOSSSHSSSSSHHHSHSHOSHHHOSSSHOSOOSSSOOOSHOHSSSHOSESOOEESCOOCE 
| ¢ WRITE DATA SET ON DISC 
j COSSSSSSSOSHSSOES SESH S SEH OS OSHS OOHSOOHHE SOOO EOOHESEOHOOOOEOEOSEOOEE 
273 MYFCB (4) =4 

274 IF (ISET.EQ.2) MYFCB (4)= 8 


275 C #S2#S2#S7% NOTE BINRY WRITES ARE a ae er ENTIR SECTR FOR SPEED #82 
WRITE(13) ((DATA (J,K),J=1,60), K=1 


aqaan 


nen nnn Ot ee ene en ae 


Cc OR eit ata TO STATOS 
ond PLEASE NOTE THE VARAIBLE SUBSCRIP PROBLEM INPUT &0 


Cc IE., VDM FORTRAN WILL NOT ACCEPT IMPLICIT DO pos IN THE READ 
WRITE STATMENT WHEN THE FIRST ELEMENT IS VAR 
Gosenennccedoucsessbsecsacososesausetersstaecteressencecrteretesébee 
PEE CLISTLD , | 11) NCODE(NSET) ,NCODE (NSET+1) ,NCODE(NSET+2) , 


SBBRIRRIRE SIG 


I A! ie BE iy 


88 


WWWWWWWWWW 
eee ee ee 
WOBNAUNAWN-—O© 
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11 FORMAT('1* ,4A2) 
WRITE(LISTLU, 4) (IPLOT (1) ,1=*48,60) 
WRITE(LISTLU,12) (NOTE(I), I=1,35) 


12 FORMAT (* *,35A2/) 


WRITE(LISTLU, 13) 


13 Poane seas T4, 4 Bite Sar, tay ris, *EZRO’, 


T28, *TRO *THHO 
‘URITECLISTLU, 14) PARAM(INDEX) , PARAM (INDEX* 1) PARAM ( INDEX+2) , 
1 PARAM (INDEX+3) , PARAM ( INDEX+4) , PARAM ( INDEX+5) 


14 FORMAT (1X,@6PF7.3, Ti@,@PFS.2, T16,@PF7.2, 124,1PE11.4, 


1 T36,1PE11.4, ‘T48, 1PE11.4//) 


WRITE(LISTLU, 15) 
1S FORMAT(1X, TS,’ML', 1T13,°’EH’, T22,’PH’, T3@,’*MCMOH’) 


WRITEC(LISTLU,16) ((DATA (J,K), K=1,4), 


? FORMAT(1X, 1T2,QPF7.3, 110,QPF7. 2, re, “OPFe i 727, 1PEL1.4) 
Crate sae We erieial a ai sen hacey ol iO wan auare tiem sia scaialerate GET SECOND IF NOT IN i.7 


QAAQAAAAAANANAAAQ 


elo le lolo i © i © i © mo] 


IFCISET.EQ.1) GO TO 2 


PRELIM PROCESS FINISHED : 


SROSSSSSSSSSSSSHSHKE SSE SSHPSCSOOSSHSSHSHSEHOSHSOSSOHASTSSSOSSHOOHOOSEECE 


DO LINEAR NEAREST NEIGHBR 
INTERPOLATION IF DEFINED 
COOCCEOCOSEOOCEEEEEEEEOOUO EEO OEEOSOOEOEOOCOOEEOHEOSOOHOOHOOOOOOOCE 


RESET THE "LOWER BRACKET FOUND" FLAG TO NOT FOUND CONDITION. 
INITILIZE THE NEXT POTENTIAL INTERPOLATABLE POINTS INDEX PARAM 
SET THE INDEX POINTER TO THE SECONDARY DATA SETS START 
Now 0 SAME FOR THE PRIMARY DATA SET 


a ee oe we ee = 


COMPARE THE Ay DATA SETS INDEPENDENT VARIABLE TO THE CURRENT 
SECONDARY SETS INDEPENDENT VARIABLES VALUE. 


IF A LOWER BRACKET HAS BEEN SET, IT IS THE CLOSEST POINT 
AVAILABLE IN THE FURNISHED DATA SET. SINCE THE CLOSEST UPPER 


hd 


SHESSESSBAR AAAS SBR RR BRISA SSSR aaE 


7 
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Cc UPPER BOUND, HAS ALSO BEEN FOUND HENCE WE'VE FOUND THE MINIMUM 
Cc INTERPOLABLE INTERVAL. GO INTERPOLATE. 
6 IF (IBRAKT.EQ.1) GO TO 19 
C 
C NO LOWER BOUND HAS YET BEEN SET, HOWEVER IF THE PRIMARY AND 
: =e SETS ARE EQUAL WE HAVE DEGENERATELY INTERPOLATED, SO 
x IF(DATA (J,5).BQ.DATA (K,3)) GOTO 20 
C 
C NO INTERPOLATIONS WERE POSSIBLE OVER ENTIRE ae IF NEXT TRUE 
- IF((J.EQ. IPLOT(39)) . AND. ((L-1) .EQ.@)) GO TO 100 
C SOME _INTERPOLATIONS WERE POSSIBLE, BUT NEAR THE END THEY WERE NOT 
Cc HENCE SAVE THOSE WHICH 3 
IF(J.EQ. IPLOT(39)) GO TO 2 
C ELSE TRY PRIMARY SETS NEXT 1 NDEPENDENT VALUE, NO BRACKET HAVING 
Cc MATERIALIZED. 
J=J+1 
Go TO 17 
Cc 
Cc 
& 
. IN THIS LOOP WE’RE LOOKING FOR AN UPPER BOUND FOR INTERPOLATION 
C EXIT_IF LAST SECONDARY’S INDEPENDENT IS LESS THAN THE PRIMARY’S 
c FIRST INDEPENDENT POINT (IE., NON COINCIDENT DOMAINS) 
¢ 18 IF((K.EQ.IPLOT (4@)).AND. ((L-1).EQ.0)) GO TO 1001 
C IF THERE ARE NO MORE ah ee POSSIBLE, BUT SOME HAD OCCURRED 
Cc DUMP THE GOOD STUFF T 
e IF (K.EQ. IPLOT (49) ) co TO. 3° 
C 
C BUMP ny. te DATA SET’S INDEX 4 NEXT TEST,AND SET 1 geiciaalaaiaas 
Cc THAT A LOWER BOUND WRT THE PRIMARY POINT IN QUESTION 
c HAS BEENFOUND. 
K=K+1 
IBRAKT=1 
GO TO 17 
Cc 
Cc 
Cities Seca cent ena anomie naa Soe ea Rone eee modcennt 


a a a a ee = 


19 DATA (L,2)= (DATA (J,5)-DATA ((K-1),3)) © (DATA (K,4)- 
Pog ae / (DATA (K,3)-DATA ((K-1),3)) + DATA ((K-1),4) 


c 

C 

C “DEGENERATELLY INTERPOLATE — 
‘ 20° DATA (L,2)=DaTA (K, 4) 
C 


.: 87. 


] 
7 
4 
i 
| 
| 


4 PAGE 8 08/09/77 VORTEX CTUF 824 HOURS 
a, 
407 (-------------—-----~-------------------------- +--+ ----- +--+ - =~ - + 
: 408 C CALCULATE THE DIFFERENTIAL DEPENDENT (IE. SEC-PRI) VALUE 
469 Cc ON THE INTERPOLATION, THEN NORMALIZE WRT FACTOR PARAM (2) 
: 419 C WHICH FOR THE CASE AT HAND IS THE TOTAL ORGANIC MATTER CONTENT 
i 412 21 DATA (L,1)=(DATA (J,6)-DATA (L,2))/PARAM (2) 
41 
ise 
t 418 c STORE ALL DATA CORRESPONDING TO A DIFFERENTIAL MEASUREMENT 
41 ee en ee eee 
418 Cc 
419 € PLACE PRIMARY INDEPENDENT INTO THE RESULT MATRIX 
i 420 : DATA (L,5)=DATA (J,5) 
4 
422 C PLACE THE PRIMARY DEPENDENT INTO SAME ARRAY 
= és DATA (L,6)=DATA (J,6) 
: 425 C IF INTERPOLATIONS HAVE BEEN ATTEMPTED OVER A 
426 (>. PRIMARY INDEPENDENT VARIBLES, THEN WE’ RE DONE. GO DUMP. 
pel ‘ IF(J.EQ.IPLOT (39)) GO TO 22 
4 429 c ELSE BUMP THE INDEX TO THE NEXT PRIMARY INDEPENDENT VALU. 
+ =Jj+ 
431 c 
432 Cc AND INCREMENT THE COUNTER HOLDING THE # OF VALID PEEP CLAT ORS 
<= Cc Be HAVE BEEN POSSIBLE TO THE NEXT POTENTIAL VALUE 
‘ 4 =L+ 
435 GO TO 17 
436 Cc 
! 437 c 
. 438 c : 
439 (----------+----~---~---~---~------- ++ --- +--+ -- +--+ +--+ 
440 . NOTIFY OPERATOR NO INTERPOLATION POSSIBLE 
: 44) See nee ee ew ee eee een we ee ee ewe eee wee sees eeeeee eee eee ees eeeesoorooeoonau 
; 442 1001 WRITE(LISTLY, 1000) 
¥ 443 10@@ FORMAT(* NO INTERP POSIBLE’) 
444 GO TO 26 
445 c 
446 Cc 
poli Meaeebasenciieb ruavanecceacuvannvsahdbaentbeeambesebakacwiaesitin 
449 c DUMP FINAL CALC(DATA 
45@ Cc 2 B PLOTD) TO DISC 
451 CHSSHSSSOSSHOSOSHESOSHV SESE ESOS HSSSOEOSHSHSSSOOHOSSOSHHSSHOOSOHOOOOOOSLOEOE 
452 22 MYFCB (4)= 12 
453 WRITE(13) (CDATA (J,K), J=1,60), Ke1,2) 
454 MYFCB (4)= 14 
455 WRITE(13) ((DATA (J,K), J=1,60), K=5,6) 
456 Cc 
4s? ] 
458 C 
poe Nsntnghdedekeuaibacinewiiaetaeeinaieeudienscietucdsaebsieainds 
461 Cc SET UP PLOT REQ SECTO 
462 Cc (WE ARE PREPARING TO USE BIPLOT) 
462 COSSSOHSSSSOHSSSSEHESESSSESOEHEOOHESHOSSOSSHESHESSEEOPOUVHSOHHESCHOOOOHOCEE 
464 G 


SESsRE8 


AMNAAAAAAANS 
et et ee et et eet et et ee 
ODAIDNALWN—-O© 
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lol oleh) 


FIRST X 
NOTE L IS THE INDEX OF THE INTEROP- 
ATABLE DATA POINTS ARRAY SEE 
SOURCE REGION 300 , LININTERPOL. 
IPLOT (9)= 14 
IPLOT (10)= L 
IPLOT (11)=+] 
Cc THEN Y1 
IPLOT (12)=12 
IPLOT (13) =L 
IPLOT (14)= +4 
c SET END Y AXI FLAGS 
IPLOT (15) =0 
IPLOT (16)=0 
IPLOT (17)=0 
Cc WE WANT LEDG AND BX&LN 
IPLOT (1)=1 
IPLOT (3)=+1 


ee a ew rw ww on ww oe we ew oe eo ee 


Cc 
C DUMP PLOT REQ SEC T DISC 
. (REMEMBER WE POINT TO 3) 
MYFCB (4)= 3 
WRITE(13) RPLOT, IPLOT, LABLX, LABLY 


SHOW HIM PLOT REQ 
FOR FUTURE REFRNCE 
WRITE(LISTLU , 4334) 
4334 FORMAT(’ {PLOT REQ SECT’) 
WRITE(LISTLU ,8818) RPLOT 
$818 FORMAT (IPSE1@.3) 
WRITE(LISTLU,8819) (IPLOT(I), I=1,47) 
8819 FORMAT (1X, 12(16, 1X)) 
WRITECLISTLU,12) (IPLOT (1), 1=48,60) 
WRITE(LISTLU, 12) LABLX 
WRITE(LISTLU, 12) LABLY 


a a errr 


DISC FORMAT POST CRUNCH....... THEN SHOW HOW HIS DATA 
DISC AREA WAS WRITN 
WRITE(LISTLU ,6543) LUN, (MYFCB (I), I=8, 10) 
6S43 FORMAT (///* BINRY DATASECTRS A. 44 ; LUNs! ,I3, 3X, ‘FIL=', 3A2/ 
1=SYSTM' 3X, '2=SYSTM’ ,3X, °3=PLOTR 
, 4=ML PRI, 3X," S=EH PR* “ax° *6=PH PR’, “3x, *7=PMCMOH’ / 
* SeML SEC' 3X. ‘9=EH SE’ ,3X,'1@=PH SE’, 3X, ’11=SMCMOH’/ 
* 12=MCMPOC’ ,3X, * 13=SMCMOH® , 3X,'14=PH'PRi (INTERPOLATABLE) */ 
* 15=PMCMOH (F(PH PRI))") 


Cc 
C 
COSSOSOSSOSSSOSSSSSESSSOOHSOSOOOSEDESOOSHEHOSHLOSHHSSEOOOSOOHOOOHOEOE 
Cc 
C 


DUMP DAT 2B PLOTD 
TO STATOS 


89, | 
" 
1 ] 
i} | 
} 
ah PAGE 10 08/09/77 VORTEX CTUF $24 HOURS ] 
5 | §23 CHSSSHEHSSSSSHSSHSSSSSHSHSHSOSSOHSHOSSOSSSSSSHSOSOSSSHHOSSOSSSOSOSOOOOOOE ] 
ji $24 WRITE(LISTLU, 23) 
525 23. ~FORMAT('1 °) 
328 ‘ WRITE(LISTLU, 4) (IPLOT (I), I=#48,60) 
528 WRITE(LISTLU, 24) | 
$29 24 FORMAT(1X/' ',T12,°H ACTIVITY'/ TS,°PH’, T14,°MCMPOC’, | 
$30 ‘ 1 T26,’PMCMOH', T38, *SMCMOH') | 
| 
‘ 532 DO 26 J=1,L 
533 WRITE(LISTLU, 25) DATA (J,5), DATA (J,1), DATA (J,6), DATA (J,2) 
34 3 FORMAT(1X,@PF8.4, T12,1PE10.3, T24,iPE1@.3, 136, 1PE10.3) 
\ 536 OG 
' §37 Cc 
} 538 Cc 
i $39 Cc 
§ 540 26 CONTINUE 
i $41 C KICK OUT LAST PAGE 
$42 WRITE (LISTLU, 1002) 
; $43 1002 FORMAT(' ®®*CRVDIF BYE®*e’) 
i 544 WRITE (LISTLU , 23) 
545 oN 
c 346 Cc 
; 246 Nu Aieiasele allen daw eaeeeecenensdeeeebadieenheceasanaeannansennenees 
549 c CLOSE PLOT FILE OBNE UPDATE 
5580 CHSSSSSOSSSSSOSSSHSHSHSEHHOHSSSHSHSSSHSOSSSHSSOSESESHSOHOSEOSSEOOSOOCOSE 
j $51 CALL VSCLOS (13,0) 
F| 552 c RETURN TO ROOT SEGMENT 
) 553 c OMTITR (IF AS OVLAY) 
554 RETURN 
| 555 END 
556 
' 557 
558 "Oro 


at ed 


a a 


PAGE 


ESSLSSRSRREBBLSBSSBRROHLISSIBALGH=Sowr)9MHMEuNnH 
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CHHHHRREREARERHHAHHHAEAMRR AA HHAAER RAHA AHREMHRHRERE RNR RREERRRR HERERO HHH 
CHHEBHHHRARAHAHHRE EMAAR EMR RAR RA RRHARHH EKA 
EEE ddd dididdddd dud dl ddudiddediddeddadad 


NLNOM ~ VI77 WITH DISC WRITE HUIZENGA 


Cc 
CHHAHRHEHHAHHEHAIAH HANAHAN MHA AA AAR RRR AAR 
ee kd dd ddd dadidddadudidididididuddududiddadldudidudadadadad 


ANAQAAAANAANANAANAAAAAN 


OH CONSUMED BY OM 
TOTAL # OF SITES C2=TMCMPOC 


NLNOM - A NONLINEAR LEAST-SQUARES FIT TO OM oe DATA 
WITH A FOUR PARAMETER MODIFIED HENDERSON-HASSELBACH 
RELATIONSHIP. 

DEPENDENT VARIABLE X IS THE MCMPOC (OH CONSUMED BY OM) 


X=C2/(10°* ((PK-PH(1))/N)+1)-Cl 
DERIVED FROM RELATIONSHIP 
PH=PK+NLOG!0(ALPHA/ 1~ALPHA) 
WHERE ALPHA IS DEGREE OF DISSOCIATION 


IN THE FORM APPLIED TO THIS DATA 
PH=PK+N®LOG( (X+C1)7(C2-(X+C1))) 


X=MCMPOC 
C1l=SHIFT OF CURVE 


CHHHHHHHHHHHHHRHHAHH RARER REMHAHAH RRR 
eo ner ede ihe a ST eee 


Cc 


agaaAaN 


DISC WRITE TO LUN=23 SECT 4&8 FILE=IBINRY 


TITLE NLNOM 

DIMENSION SUMSQ (2) ,PH(6@) ,X(6@), IFCBUD(13) 

DIMENSION IDUMMY(120), AMODEL(66), ADIFF (60) 

bide a (3)/7* °/,(IFCBUD (1), 1=8,10)/’ IB’,’ IN’, *RY’7 


GET DATA AS STORED BY TCRVDI ON FILE IBINRY 


CALL Sate IFCBUD , IREW) 


WRITE (1 , 160 

FORMAT (’ 1 THIS IS NLNOM’) 
WRITE(1, 102 

FORMAT (’ INPUT LISTLU, IOLUN’) 
WRITE(1, 106) 

FORMAT (’ +### +###") 


READ(1,168) LISTLU, IOLUN 
FORMAT (1X,14,1X,I 
WRITE(IOLUN , 160) 
WRITE(LISTLU , 100) 
WRITECIOLUN, 109) 


AD-A0S5 157 RHODE ISLAND UNIV KINGSTON GRADUATE SCHOOL OF OCEANO=-ETC F/G 
SUPPLEMENTAL DATA REPORT ON INVESTIGATIONS OF MARINE DISSOLVED wnerewu) 
MAY 78 D L HUIZENGA N00014-76-C=-0226 : 
UNCLASSIFIED URT/GSO-REF=-7A-1 we, 


Ce 2 END 
AD ; DATE 
ADBS It FILMED 
7 ~78 
iE] Doc 


91. 
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59 109 = FORMAT(* Jed gr ene ine VALUES: C1 C2 PK N’) 
WRITE(IOLUN, 1 

61 11@ FORMAT(’ +##. tone +hH HEH FHM HHH HH OOOR!) 

62 READ(IOLUN, 112) C1, C2, C3, C4 

o 112 FORMAT(IX,F8. ee ms F8.4; 1X, Fs. 4,1X,F8.4) 

65 

66 


lanes he SAR AI ae ee 


WRITE CIOLUN, 
114 FORMAT(’ WHAT PH RANGE?’ ) 
] WRITECTOLUN, 116) 
4 67 116 FORMAT (1X, WH, HHH HH HHH’) 
i 68 READ(IOLUN, 118) PHMIN, PHMAX 
; 69 118 FORMAT(1X,F7.3,1X,F7.3) 
i 70 WRITECIOLUN, 120) 

71 126 FORMAT(* CONVERGENCE CRITERION’) 


72 WRITE (1OLUN, 122) 
; 73 122 FORMAT (1X, +####8### IN E FORMAT °) 
; 74 READ (IOLUN, 124) CONV 
75 124 FORMAT (1X, E8. 5) 
, 76 WRITE (IOLUN , 121) 
77 121 FORMAT(? WRITE MODEL CURVE TO DISC? 1=¥ES'/1X,+#") 
78 READ(IOLUN , 123) IDISC 
79 123 FORMAT(1X, i2) 
80 SUMSQO=- 1. 
81 K= 1 
82 WRITE(LISTLU, 126) CONV, PHMIN, PHMAK 
83 126 FORMAT COITERAT IVE. PROCEDURE CONV=’,1PE13.6,’ PH RANGE ’, 
84 +OPF7.3,' - ',OPF 
i 85 WRITECLISTLU, 128). K,C1.C2,C3,C4 
4 86 128 FORMAT(’@’,16,4(1X, iPEI3.6)) 
87 158 C=@.434294 
88 FIF1=0. 
89 F1F2=0 
90 F1F3=0 
91 F1F4=0 
i 92 F2F2=0 
| 93 F2F3=0 
i 94 F2F4=@ 
95 F3F3=0 
96 F3F4=9 
97 F4F4=0 
98 CFO=6 
99 CF1=0 
100 CF2=0 
101 CF3=0 
| 102 CF4= 
103 PK=C3 
{ 1 AN=C4 
185 SUMSQ (2)s0. 
| 106 SUMSQ (1)=0. 
| i 
| 169 C NOW CALCULATE THE DERIVATIVES OF THE FUNCTION WITH RESPECT | 
110 C TO EACH OF THE CONSTANTS TO BE DETERMINRD. 
111 € CALCULATE THE ‘SUMS OF DERIVATIVE PRODUCTS. 


113 DO 200 I=1,60 

114 IF(PH(I) .LT.PHMIN) GO TO 200 
115 IF (PH(1) .GT.PHMAX) GO TO 300 
116 VAR=EXP ( (PK-PH(1)) /(C®AN) ) 


88/09/77 VORTEX CTUF 


Fie-1]. 

F2=1./(VAR¢1.) 
F3=-C2°VAR/ (C*AN® (VAR+ 1.) @@2) 
F4e-C2¢ (PK~PH(1) )®VAR® (~ 1.0) 7(C® (AN® (VAR* 1) 902) 
FeC2/(VAR+1)-C1 
FIFJ*FIF1+F1°F ] 
FLF2*F1F2+F1°F2 

F2F 1 =F 1F2 

FLIF3=F (F3+F1¢F3 

FSF 1=FIF3 

FIF4=FIFA+F 1°F4 

F4F1sF1F4 

F2F2=F2F2+F2°F2 
F2F32F2F3+F2°F3 

F3F2=F2F3 


F2F4=F2F4+F2¢F4 
F4F2=F2F4 
F3F3*F3F3+F3°F3 
F3F4="F3F4+F3¢F4 
F4F3=F3F4 
F4F4sF4F4+F4°F4 
CFO=CFO+ (X (1) -F) ©%2 
CF L=CF 1+ (X(1)-F) ®F 1 
CF2=CF2+(X(1)-F) @F2 
CF3=CF3+ (X(1)-F) ¢F3 
CF42CF4+(X(1)-F) @F4 


83 


-LT.0.) GO TO 3 
.GE.CF@) GO TO” 301 


NSQO-CFO) /CFO) .LT.CONV) GO TO 1600 
WRITE(LISTLU,302) K,C1,C2,C3,C4,CFO 
FORMAT (1X, 16,5(1X, 1PE13.6)) 


e= 


ALGRBRA TYPE CALCULATIO 
TRIX OF THE GENERAL BOUATIONS se | NO REARRANGSMENTS 
ARE MATRICES WITH RIGHT HAND SIDE OF THE SET 
OF + ag SUBSTITUTED FOR THE 
COLUMN N. 


T ALL THIS DOES IS TO SOLVE THE SET OF FOUR EQUATIONS 
BY DOING CALCULATIONS ON 4X4 MATRICES. 


BINAAAANANAANAANN 


10)» De FIF1¢ (F2F2¢ (F3F3¢F4F4~-F4F3*F3F4) -F2F3¢ (F3F2°F4F4-F4F2°F3F4) + 
+F2F4® (FSF2°F4F3-F4F2°F3F3) ) 
+<Pif2e (PoP 16 (PaR9¢F4P4-FAPS*F3F4) -F2P9* (FOP 1° F4P4-FAF 1*F3F4) 
++F2F4° (FSF | °F 4F3-F4F 1°F3F3) > 
++F 1F3° (F2F | © (F3F2¢F4F4-F4F2°F3F4) -F2F2° (FSF | °F4F4-F4F } °F SF 4) 
++F2F4° (FSF 1¢F4F2-F4F 1°F3F2) ) 
LattsiuaGs: CC ee 
++F2F3® (FSF 1°F4F2-F4F 1 ®F3F2) ) 
Dis (CFs (RoF2e (Fghgeh4ha-R4P 3° F3P4)-P2F3* (FOP2¢F4P4-PAF2¢F 34) + 
+F2F 4° (P3F2°F4F3-F4F2°F3F3) ) 
o-F1F2¢ (CF2¢ (F3F3°F4P4-F4F3°F3F4) ~F2F3°® (CF3¢F4F4-CF4°F3F4) 


eT ee 
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od 
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330 


LAAQQONMONO 


Es 


aia 


++F2F4¢ (CF3°F4F3-CF4°F3F3) ) 

++F 1F3° (CF2¢ (F3F2°F4F4-F4F2¢F3F4) ~F2F2¢ (CF3°F4F4-CF4°F3F4) 
++F2F4°¢ (CF3°F4F2-CF4°F3F2) ) 
+-F 1F4¢ (CF2¢ (F3F2°F4F3-F4F2°F3F3) -F2F2¢ (CF3®F4F3-CF4°F3F3) 
++F2QF3® (CF3®¢F4F2-CF4°F3F2) ))/D 

eet itatcatetatia eee (CF3°F4F4-CF4°F3F4) + 
+F2F4° (CF3°F4F3-CF4°F3F3) ) 

+-CF 1° (F2F 1° (F3F3¢F4F4- Seana ot a) -FSES* CFSE IOFAPA-FAP 1SR SFA) 
++F2F4° (F3F 1 °F4F3-F4F 1 °F3F3) 
+ SPIPS* (FSF 1* (CE3*F4E4-CR4® PSF 4) ~CF2° (FSF 1°F4F4-F4F 1°P3P4) 
++FQF4® (F3F 1*CF4-F4F 1 *CF3) ) 
+-FLF4® (FOF 1° (CF3®F4F3-CF4°F3F3) -CF2° (F3F 1 °F 4F3-F4F | °F3F3) 
++F2F3° (F3F 1®CF4-F4F 1®CF3) >) /D 

D3s (FIF1® (F2F2* (CF3¢F4F4-CF4°F3F4) -CF2° (F3F2°F4F4-F4F2°F3F4) + 
+F2F4°¢ (F3F2¢CF4~F4F2°CF3) ) 
+-F LF2® (F2F 1 (CF3¢F4F4-CF4°F3F 4) -~CF2° (FSF ] OFAF4-F4F | °F3F4) 
++F2F4¢ (FSF 1¢CF4-F4F 1°CF3) ) 

++CF1® (F2F 1° (F3F2@F4F4-F4F2°F3F4) -F2F2° (FSF 1 °F AF4~-F4F 1 °F SF 4) 
++F2F4° (FOF 1®F4F2-F4F 1*¢F3F2) ) 

+-F 1 F4® (F2F 1% (F3F2°CF4- ej ee 
++CF2¢ (F3F 1°F4F2-FA4F 1°F3F2)))/ 

Das (FIFI (Foro (FSF 3° Crd Fak 3° CF3) - F2F3° (F3F2°CF4-F4F2°CF3) + 
+CF2® (F3F2°F4F3-F4F2°F3F3) ) 

+-F 1F2® (F2F 1® (F3F3®CF4~F4F3°CF3) -F2F3® (F3F 1 ®CF4~F4F | °CF3) 
++CF2° (F3F 1®F4F3-F4F 1°F3F3) ) 

++F 1 F3° (FOF 1® (F3F2@CF4~F4F2°CF3) -F2F 2° (FSF 1®CF4-F4F 1°CF3) 
++CF2° (F3F 1¢F4F2-F4F 1°F3F2) ) 

+-CF )® (F2F 1° (F3F2°F4F3-F4F2¢F3F3) -F2F2¢ (FSF 1°F4F3-F4F 1°F3F3) 
pe sith die id i die 

1 


KeX+] 


NOW DO CALCULATION OF SUM OF Ry WITH BOTH 
THE FULL CORRECTION TO THE FOUR PAR 
DETERMINED ABOVE AND WITH ve THE: CORRECTION. 


DO 600 [*1,2 

Al#FLOAT (I) /Y 
CE1L*Ci+DI®Al/2. 
CE2=C2+D2"Al/2. 
CE3*C3+D3"Al/2. 
CE42C4+D4eAl/2. 

DO SOO J=1,60 
IF(PH(J).LT.PHMIN) GO TO SQ 
IF (PH(J) .GT.PHMAX) GO TO 600 
VAR#EXP ( (CE3-PH(J)) /(C®CE4) ) 
XVeCE2/(VAR+1) - CEL 
DIFF2X(J)=-XV 

SUMSQ (1) =SUMSQ( 1) «DIFFe®¢2 
CONT INVE 

CONTINUE 


CALCULATE THE AMOUNT OF THE CORRECTION TO BE MADE WHICH WILL 
YIELD THE LOWEST SUM OF SQUARES. 


GN Ch th 


94, 


68/03/77 VORTEX CTUF 831 HOURS 
% SOLUTION TO QUADRATIC EQUATION GIVES AMOUNT TO APPLY. 
Cc 
. VMIN=0.5/¥+0.25/Y¥* (CFO-SUMSQ (2))/(SUMSQ (2)-2.°SUMSQ (1)+CF®) 
Cc IF CALCULATED VMIN IS LESS THAN ® GO BACK AND DO THE 
c CALCULATIONS AGAIN ONLY HALVE THE INTERVAL OVER WHICH 
: THE CALCULATIONS ARE MADE. 

c 
IF(VMIN.LT.®.) GO TO 700 
SUMSQO=CFO 
C1#C1+VMIN®D1 
C2=C2+VMIN®D2 
C32C3+VMIN®D3 
C4=C4+VM INO DS 
GO TO 1S@ 
re YeY®2. 
£ ONLY ALLOW THIS SHORTENING OF THE INTERVAL TO GO SO FAR 


i a GO TO 1800 


GO TO 460 
1000 WRITE(LISTLU, 1001) K 
1601 FORMAT('@CONVERGENCE CRITERION MET NUMBER OF ITERATIONS=’ , 14) 
Cannan ear eer eet ter ae anne teen con cabernet ree eneeeee 
Cc E TO DISC THE MODEL CURVE AND RESIDUALS 
Cc LUN=23 FILE=IBINRY 
C SECT 4 pe ® pens ee 
C SECT 8 RESIDU 
CHHMADEDNEANSUSANSHNDRAYSASERNERAOAHHHRAHANUMRHERONHHNOROEERHEOD 
DISC.NE.1) GO TO 1200 


is) 

1069 DO 1100 I=*!,NUMBER 

VAR=EXP ( (C3~PH(1))/(C®C4) ) 

AMODEL (1) =C2/(VAR+1)-Ct 

ADIFF (1) =X (1)-AMODEL (1) 
1108 CONTINUE 

IFCBUD (4) =4 

WRITE(13) AMODEL 

IFCBUD (4) =8 

WRITE(13) ADIFF 
12@@ CALL V8CLOS(13, IREW) 
OTT Addl ddd didodudududodudadidodododododudodedodidudidodedadidocededededsdubsdudedpdidodudodudodidudududadududodedua 
aca ie kel a ihe a la ea dial 


Cc OUTPUT FINAL PARAMETERS 


waerrie aie te Cl 
MAT ¢ Cle ',F8.4) 
WRITE CLISTLU, 1003) C2 
C2 .4) 


FORMAT ( = '.F8. 
WRITE(LISTLU, 1006) CFO 
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1006 FORMAT(’ SUM OF SQUARES «= ',1PE13.6) 
WRITE(LISTLU, 1008) NUMBER 

1088 FORMAT (1X,'# OF POINTS IS ',14) 

WRITE(LISTLU, 1607) 

WRITECIOLUN, 1007) 

FORMAT (* PePeNLNOM BYES? *) 


1007 
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Copy of cards used to run the SAS 76 NLIN program for non-linear 


least squares model of combined titration curves. 
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